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ABSTRACT 

We present a new assessment of the ability of Infrared Dark Clouds (IRDCs) to form massive stars 
and clusters. This is done by comparison with an empirical mass-size threshold for massive star 
formation (MSF). We establish m(r) > 870 M Q (r/pc) 133 as a novel approximate MSF limit, based 
on clouds with and without MSF. Many IRDCs, if not most, fall short of this threshold. Without 
significant evolution, such clouds are unlikely MSF candidates. This provides a first quantitative 
assessment of the small number of IRDCs evolving towards MSF. IRDCs below this limit might still 
form stars and clusters of up to intermediate mass, though (like, e.g., the Ophiuchus and Perseus 
Molecular Clouds). Nevertheless, a major fraction of the mass contained in IRDCs might reside in 
few 10 2 clouds sustaining MSF. 

Subject headings: ISM: clouds; methods: data analysis; stars: formation 
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1. INTRODUCTION 

About a decade ago, Galactic Plane surveys re- 
vealed large numbers of Infrared Dark Clouds (IRDCs, 
lEgan et alj|1998t iPerault et all 11996ft . These are iden- 
tified as dark patches against the diffuse Galactic mid- 
infrared background. First studies of very opaque IRDCs 
suggested that these have very high densities, column 
densities, and masses ( nfH^] > 10 5 cm -3 , A^ty > 
10 23 cm^ 2 , m > 10 3 M^: ICarev et al.lH99ll . Since they 
are dark, they are likely to be in an early evolution- 
ary phase. Embed ded in IRDCs are "cores" of a few 
dozen solar masses (jCarev et all 120001 ) . It has therefore 
been suggested that many IRDCs are the long-sought 
examples of clouds just at the onset of the formation 
of massive stars and (proto-)clusters. This notion was 
corroborated by observati ons of young massive s t ars in 
a few individual I RDCs dRathborne et all [2001 [20071 : 
iPillai et all [2001 iBeuther fc Steinackerl 12007^ Such 
views also f orm the framework of schemes for IRDC evo- 
lution (e.g.. lRathborne et al 1 120061 [ Rygl et al.ll2010D an d 
reviews (e.g.. iMenten et alj|2005l IBeuther et all 120071) . 
IRDC samples are usually compared to regions of mas- 
sive star formation (MSF), such as Orion and M17 (e.g., 
iRagan et alll^OOl ). 

This pict ure cannot be complete , though. The above 
studies fand lPeretto fc Fuilerll2009l ) acknowledge that re- 
gions forming low and intermediate mass stars can also 
appear as shadows in images at mid-infrared wavelength 
(jAbergel et al.1 119961 ). Such IRDCs will not form mas- 
sive stars. Unfortunately, the number of IRDCs evolving 
towards MSF is presently not known. Fractions up to 
100% have been considered in the past ( Section 14. 3|) . 

In this letter, we thus use a novel criterion to provide 
the first conclusive quantitative demonstration that only 
few IRDCs are headed towards MSF. This aids identify- 
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ing pre-MSF IRDCs as targets for ALMA and Herschel. 
As a bonus, the MSF threshold identified below — the first 
observational limit of this kind — informs the ory. 

In papers I and II (|Kauffmann et al.H2010allbh . we show 
that solar neighborhood clouds devoid of MSF (specifi- 
cally: Perseus, Ophiuchus, Taurus, and Pipe Nebula) 
generally obey 



m(r) < 870 M (r/pc) 



1.33 



(1) 



IRDCs submitting to Eq. ([1) would resemble, e.g., Ophi- 
uchus and Perseus, but not Orion (which violates Eq. 
[]}. Figure [T] illustrates why clouds bound for MSF must 
exceed Eq. (TTJ). Since star formation necessitates an ap- 
propriate mass reservoir, MSF requires that a large mass 
is concentrated in a relatively small volume. Based on 
more detailed theoretical considerations, Section l4~T1 puts 
quantitative limits on this intuitively evident reasoning. 
As seen in Fig. [IJ the masses in this MSF region are well 
above the mass-size range bound by Eq. ([T]). Observa- 
tions of MSF clouds confirm Eq. (fT]) as a true MSF limit 
(Section 13. ip . This suggests to use Eq. ([1} to roughly 
separate IRDCs with (future) MSF from those without. 

This letter is organized as follows. Based on data 
from Section [2j Section 13.11 confirms (using known MSF 
clouds) that Eq. ([1]) approximates an MSF limit. Many 
well-studied IRDCs (25%-50%) fall short of this thresh- 
old (Section 13. 2\ . Less certain data for complete IRDC 
samples suggests that most IRDCs obey Eq. ([T]), and will 
thus not form massive stars (Section I3.3p . Still, most of 
the mass contained by IRDCs might be in clouds forming 
massive stars (i.e., those violating Eq. [T]). 

2. METHOD & DATA 

2.1. Sample 

Data for solar neighbourhood clouds not forming mas- 
sive stars (here: Taurus, Perseus, Ophiuchus, Pipe Neb- 
ula) are taken from paper II (and references therein). 
We rely on bo l omete r surveys to characterize MSF sites: 
IBeuther et"all (|2002[ ) study FIR color-selected MSF can- 
didates with CS-detect e d den se gas but no radio con- 
tinuum; iMucller et al.l (|2002| ) map water masers em- 
bedded in CS clumps of high bolometric luminosity (> 
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Fig. 1. — Plausible th eoret ical mass-size limits for massive star 
formation (MSF, Section l4.ll green and blue shading) in relation to 
mass-size laws (e.g., m[r] = mo • r°) observed for non-MSF clouds 
(Eq.[T] Fig. [2] yellow shading). At small radii, MSF clouds (high- 
lighted in red) must contain fragments bound by the theoretical 
MSF-limits. Depending on the interplay of slope, b, and intercept, 
mo, such clouds must also at radii > 0.1 pc be more massive than 
fragments in non-MSF clouds. 

10 3 L Q ); iHill et all (|2005h explore meth anol masers and 
ultra-compact Hn regions; iMotte et al.l ((200 7) study the 
nearby Cygnus-X MSF site (we use their 'clumps'). To 
exclude fragments not forming massive stars, we only use 
the 'Type 1' sources (CH3OH and/or Cr TCN emission, 
no re solved radio continuum) from the iBeuther et al.1 
(|2002f l survey, and ignore the secondar y 'mm-only ' cores 
(without masers and Hn regions) in the IHill et all ([2005) 
study. 

The 1RDC s amples were crea t ed u sing MSX and 
Spitzer i mages. IRathborne et using bolome- 

ters) and lRagan et all (|200a using 8 /Ltm e xtinction) fo- 
cus on clouds with stark 8 fxm contrast. ISimon et al.1 
(|2006| ) report 13 CO-based results for a ll IRDCs evident 
in the ir 13 CO Galactic Plane survey. iPeretto fe Fuller! 
(|2009f ) catalogue extinction properties for 11, 000 Spitzer 
8 /im IRDCs with unknown distances. 

2.2. Data Processing 

The mass-size data for solar neighborhood clouds are 
derived in paper II (using methods summarized in Sec- 
tion 2.1 and Fig. 1 of paper I). They are based on column 
density maps derived from dust emission (MAMBO and 
Bolocam) and extinction (2M ASS) data. Using a den - 
dogram method introduced by IRosolowskv et al.l (|2008D . 
starting from a set of local column density maxima, a 
given column density map is contoured with infinites- 
imal level spacing. Every contour defines the bound- 
ary of a cloud fragment. We derive the contour-enclosed 
mass and the effective radius, r = (A/it) 1 / 2 . Subsequent 
contours/fragments are usually nested. This defines re- 



lationships between cloud fragments, essentially yielding 
series of mass-size measurements. In Fig. [2] such series 
are drawn using continuous lines. 

To derive column densities from the extinction maps, 
we assume that column density and visual extinction 
are related by 7Vh 2 = 9.4 x 10 20 cm~ 2 (Ay /mag) 
(|Bohlin et al.l 119781 ). To combine dust emission 
and extinction observations, they must be calibrated 
with respect to one anoth er. In practice, we use 
IQssenkopf fe Henningl (|1994f ) dust opacities (decreased 
by a factor 1.5, to match observed opacity laws; Section 
4.2 of paper I) for emission-based masses. 

For comparisons, we must scale all masses to the col- 
umn density laws from paper II. Also, it is necessary to 
harmonize the different definitions of mass and size. The 
scaled data are shown in Fig. [5] 

Where relevant, we use dust temperatures suggested by 
the original studies. However, we substitute our choice of 
dust opacities and the aforemention ed 1.5 scaling factor . 

CO masses are directly taken from I Simon 
since their 13 CO-to-mass conversion law is in rough 
agreement with (i.e., by factors of 1.1-2.0 larger than) 
the ex tinction-calibrated ones derived by iPineda et al.1 
(|2008| ). We assume that dust emission at 1.2 mm wave- 
length and optical depth at 8 [iva wavelength are re- 
la ted by F,^ eam = 50 mJ y(ll" beam)" 1 • t 8 Mm (Eq. 4 
of IPeretto fe Fullerl 120091 ) . and derive column densities 
from these intensities (assuming dust at 15 K, and using 
the 1. 5 scaling factor). We thus increase the lRagan et al.l 
(|2009l) masses (from their case 'A') by a factor 1.47 
(to account for their choice of opacities and molecular 
weightsjE 

In many cases (IBeuther et all 120021 IHill et al.l I2005L 
IRathborne et al.ll2006l IMotte et al.ll2007t) . the size listed 
in the original publication refers to the contour at half 
peak intensity, while the mass measurement includes 
emission at much lower levels. In these cases, we as- 
sume that the sources have a near-Gaussian shape (just 
as explicitly assumed in many of the original papers) . For 
such sources, the mass contained in the half peak column 
density contour is ju st a fraction ln(2) 0.6 9 of the to- 
tal mass (Eq. A.23 of iKauffmann et al.l[200l the area at 
half peak intensity is 7t[6*fhwm/2] 2 ). Thus we reduce the 
mass to a fraction ln(2), and use half of the published 
FWHM size as the effective radius. iMueller et all (|2002l ) 
list masses for a sphere, not an aperture, and so the mass 
(taken for the smaller of their ra dii) has to be scaled up 
by a factor of order n/2 « 1.57 (Kauffm ann et aLll2008l . 
Eq. 13). If more than one distance is listed for a given 
object, we adopt the smaller one (yielding a lower limit 
to m[r]/mii m [r] derived below). 

3. ANALYSIS 

3.1. A Threshold for Massive Star Formation? 

At given radius, a cloud fragment (i.e., clump, core, 
etc.) can be compared against Eq. ([1]) by deriv- 
ing the mass ratio m(r) / my ml (r) (where mn m [r] = 
870 Mq [r/pc] 1 - 33 ), to which we refer as the 'compact- 

5 We further correct their masses by factors 4/(7r[(r/pc) ■ 
206, 265"/(d/pc)] 2 /l'.'2 2 ), where d is distance, since pixels per 
beam (as erroneous ly adopted) have to be replaced by pixels per 
clump in Eq. (5) of lRagan eFaLl 1(20091) . 
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Fig. 2. — Clouds with (panel b) and without massive star for- 
mation (MSF; panel a), compared to IRDCs (panel c). Concep- 
tually, data in panels (b) and (c) correspond to picking a single 
mass-size measurement along one of the lines shown in panel (a). 
The MSF and non-MSF clouds suggest that the indicated limiting 
pow er la w (Eq. [TJ approximates a mass-size limit for MSF (Sec- 
tion [37TJ . Onl y a f raction of the IRDCs exceed this MSF limit 
(Fig. [3] Section |3 . 2D . If a star- forming region contains more than 
one fragment (i.e., clump, core, etc.), the most compact fragment 
(i.e., with maximum mfrj/miimfr]) is highlighted by a circle. 
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Fig. 3. — Percentiles of the maximum compactness per cloud, 
max[m(r)/mi; m (r)], for various cloud samples. For a given sample, 
the ratios below which a certain fraction (e.g., 25%) of the sample 
members resides are indicated by bars. Local non-MSF clouds 
(Fig. \2\& \ ) have a com pactness < 1 (Eq. [TJ. The bars for the 
I Rathbo rne ct al. (2006) sample include (left), respectively exclude 
(right), their 'em' cores. Clearly, the iRDCs do not reside in the 
mass-size space unambiguously associated with MSF. 



ness' . "Secondary co res" (only listed by iRagan et al.l 
and iRathborne et al.l ) are suppressed by characteriz- 
ing star- forming regions (i.e., a given massive star, 
or an entire IRDC) by their most compact fragment, 
max[m(r)/mi im (r)]. 

Figure [3] gives max[ro(r)/TOi; m (r)] as derived for the 
samples examined here. This is based on the mass-size 
data presented in Fig. [2j The compactness assumes a 
range of values in every sample. This spread is captured 
by plotting several percentiles. 

As suggested by Fig. [2] we can clearly see in Fig.[3]that 
regions forming massive stars are, at given radius, more 
massive than the limiting mass, rnn m (r). In all surveys 
of MSF regions, > 75% of the clo uds have a maximum 
compactness > 1.7. One survey (jBeuther et al.1 I2002T ) 
contains a very small number of regions (~ 10%) less 
compact than required by Eq. ([1]). These regions might 
be interesting targets for follow-up studies. In general, 
though, this analysis corroborates the hypothesis that 
Eq. ([1]) approximates a threshold for MSF. 

3.2. Are IRDCs unusually Dense and Massive? 

Figure [3] provides a compact ness analysis for IRDCs . 
We separately characterize the IRathborne et al.l (|2006f ) 
sample including and excluding their 'em' cores with 
associated 8 /im sources (which are not dark). "True" 
IRDCs will have properties in between these extremes. 
Two interesting trends manifest in these m(r) / my lm (r) 
data. 

First, IRDCs have masses which are, for given size, 
comparable to those of solar neighborhood clouds not 
forming massive stars (e.g., Ophiuchus and Perseus). In 
all samples, > 25% of all clouds have a compactness < 1. 
Except for the IRathborne et al.l (|2006l ) clouds, > 75% of 
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all targets exceed Eq. (|T|) by a factor < 2. 

Second, IRDCs are less compact than regions form- 
ing massive st a rs. For example, excluding the 
IRathborne et all (|2006f ) targets, > 75% of all IRDCs are 
less compact than most (> 75%) of the MSF regions. 

In summary, the IRDCs studied here have (for given 
size) masses in betw een those of region s with and without 
MSF (where "true" IRathborne et al.l IRDCs have prop- 
erties in between the two extremes shown) . Very clearly, 
they do not reside in the mass-size space unambiguously 
associated with the formation of massive stars. 

However, before drawing final conclusions, let us 
consider some bia ses affecting our analysis. First, 
IRagan et all (12009ft de rive masses using CLUMPFIND, 
while IRathborne et al.l (f2006h use GAUSSCLUMPS. For 
the former, paper I showed explicitly that the derived 
masses are, for given radius, just < 70% of those derived 
using our dendrogram approach. For the latter, the same 
is expected, since the Gaussian fits only describe a frac- 
tion of the emission. In a given map, our characteriza- 
tion scheme from papers I and II would thus find larger 
masses. 

These biases are countered by other factors, though. 
We use the 'case A' masses (assum i ng br ight IR fore- 
grounds) provided by Ragan et all (|2009f ). Following 
IPeretto fe Fulled (|2009|), their 'case B' (fainter fore- 
grounds) appears to be more r ealistic. The masse s could 
thus be lower by a factor ~ 2 (jRagan et al.ll2009l ). Simi- 
larly, EnidaiFin (|2008f ) sugges ts 13 CO-to -mass conver- 
sion factors lower than used by Sim on et al.l (|2006f ). In 
any case, similar biases affect the data for MSF regions. 
Differences between these and IRDCs are not likely to 
only come from observational uncertainties. 

Finally, none of the IRDCs in the IRagan et al.l . 
IRathborne et al.l . and lSimon et al.l s amples are "typical" 
for the general Galac t ic pop ulation. IRagan et al.l (|2009l ) 



and IRathborne et al.l (|2006[) se lect clouds which ar e un- 
usually dark in 8 /xm images. ISimon et al.l (|2006l ) only 
characterize IRDCs which are relatively large and dark, 
and are clearly detected in 13 CO emission. All this ex- 
cludes IRDCs of low mass and density from the samples. 
Less biased IRDC samples should thus be less compact 
than derived here. 
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Fig. 4.— Like Fig.EJc), but for the lPeretto & Fuller! (2009|) sam- 
ple (projected out to various distances). 



TABLE 1 

Compact IPeretto fc FullerI (120091 ) IRDCs 



Distance" 


Number 6 


Fraction c 


Mass 6 


Mass Fraction 


kpc 




% 


1O 6 M 


% 


2 


831 


7 


2.0 


71 


1 


2218 


20 


9.8 


87 


6 


3639 


32 


23.6 


93 


8 


4778 


12 


43.2 


96 



"distance to which the sample is projected 

^number of clouds with m(r) > mi; m (r), and their total mass 

c mass and number fraction of compact clouds 



3.3. Typical Star Formation Properties in the Galaxy 

The IPeretto fc FullerI (|2009f ) catalogue lists IRDC an- 
gular sizes and column densities for the entire Galactic 
Plane covered by Spitzcr. It thus provides an ideal tool 
to derive a first idea of typical IRDC properties. Since 
they likely constitute (to our present knowledge) the typ- 
ical reservoir of Galactic star-forming gas, IRDC charac- 
teristics probably gauge the early state of Galactic star 
forming regions. 

Since no distances are known for the IPeretto fe FullerI 
IRDCs, we constrain their masses and size s assuming a 
reason able r ange of distanc e s. An alysis by ISimon et al.l 
(|2006D and Uackson et"aLl (|2008l) suggests that most 
IRDCs have distances of 2-8 kpc. Figure @] illustrates 
the derived masses and sizes, and Table Q] characterizes 
the IRDCs found to be compact (i.e., m[r]/m\i m [r] > I). 

This analysis has two interesting result s. First, by 
number, most of the lPeretto fc FullerI (|2009ft IRDCs have 
masses and sizes comparable to those of solar neighbor- 



hood clouds devoid of massive stars (i.e., they are not 
compact). This holds even when adopting the largest 
reasonable distance. Second, the compact clouds contain 
most of the mass (more accurately: most of the area- 
integrated column density) seen in these IRDCs, even 

for small IRDC dista nces. 

Unfortunately, the IPeretto fc FullerI (|2009f ) survey is 
(like most extinction studies) uncertain in the sense that 
it assumes that the diffuse Galactic emission can be re- 
liably modelled in its spatial distribution. This may not 
be true. In this spirit, the results from this section should 
be taken as an indication, not as a final result. 

4. DISCUSSION 

4.1. Mass-Size Structure of MSF Clouds 

Consider the following toy model to understand the ex- 
pected mass-size properties of MSF clouds. Stars proba- 
bly form on a timcscalc r s f slower (eg < 1) than the free- 
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fall timescak^l r s f > Tg/eg oc l/(eff (g) 1 ^ 2 )- In spher- 
ical symmetry, mass, size, and density are related by 
(g) = £ e m/(4/3 7rr 3 ), where e e < 1 takes line-of-sight 
material not associated with the sphere into account. A 
specific star formation timescale then requires that 

Further, to form a star of certain mass, M*, a mass reser- 
voir larger than M* is necessary: 

m(r) > M*/e m . (3) 

Figure [1] evaluates these limits for a star of 8M , based 
on cfhciencietQ and timescales from Spitzer observations 
of solar neighborhood clouds. Within the model, cloud 
collapse will only yield a massive s tar if initiated inside 
the bo undaries set by Eqs. J5H2]). iKrumholz fe McKei 
(|2008|) provide a similar limit, derived assuming that the 
collapsing region is heated by a cluster of low-mass stars 
(in our terminology, they use e m = 1/2). 

In order to sustain MSF, at least a few cloud fragments 
in MSF clouds must reside within the theoretical MSF 
boundaries mentioned above (Fig. [T]). The global struc- 
ture of these clouds can usually be described by power 
laws, m(r) = too ■ r b , with b < 2 (paper II). Such power 
laws imply that MSF clouds violate m(r) < m\i m (r) (Eq. 
[I]). Depending on slope (b), intercept (mo), and their 
interplay, such excesses are expected for radii 3> 0.1 pc 
(Fig. [J). This is just what we find for MSF clouds (Fig. 
13b]). 

MSF is thus only possible if a clouds' slopes are shal- 
low, intercepts are large, or both, when compared to Eq. 
(|T|). This permits a new way to quantitatively compare 
the structure of clouds with and without MSF. Pure dif- 
ferences in too imply that MSF and non-MSF clouds only 
differ in their absolute properties. Differences in slopes &, 
however, imply relative differences in the structure, such 
as deviations in the hierarchical cloud structure. 

4.2. Average State of IRDCs 

The IRDC properties mentioned in the introduction 
(n[H 2 ] > 10 5 cm" 3 , N[R 2 ] > 10 23 cm" 2 , m > 1O 3 M ) 
only seem to characterize the densest patches in very 
large and massive IRDCs. They are not well suited to 
describe IRDCs on average. 

Some IRDCs with m(r) < mn m (r) might further evolve 
and eventually undergo MSF. And particular dust prop- 
erties could, in principle, erroneously indicate m(r) < 
Tiii m (r) where the reverse is true. However, such caveats 
are not usually considered when using IRDC data to con- 
strain MSF. Thus we abstain from such considerations. 

Our study suggests that many IRDCs, if not most, are 
not related to MSF. One thus has to be prudent when 
using IRDC properties to constrain MSF initial condi- 
tions. Most studies discussing IRDCs as pre-MSF sites 
concentrated on very opaque IRDCs of large angular size. 
These clouds often violate Eq. ([1]), and many of them are 
good MSF candidates. 

6 Tff = (37r/[32 G (g)]) 1 / 2 , where G is the constant of gravity and 
{g) is the volume-averaged density 

7 The main accretion phase of a low mass star (IR-classes and I) 
typically finishes after 7 X 10 5 yr l|Evans et al,H2009t ). Typical free- 
fall timescales of their natal cores ~ IP 5 vr IIEnoch et al.l200l) then 
imply e ft 1/7. F urther, e ^ 1/2 (Kauffmann et al. 20081 Eq. 
13) and e m as 1/3 (Alvcs et al. 2007). Since massive stars might 
form faster, and the star formation efficiency is not constrained 
well, we explore 3 < r af /10 5 yr < 7 and 1/3 < e m < 1/2 in Fig.Q] 



4.3. Do most stars form in just few IRDCs? 

iRathborne et al.l (j2006f ) suggest that most of the 
Galactic star formation might come from IRDCs. The 
absence of other likely reservoirs of star-forming gas 
evinces this too. By number, most IRDCs are likely to 
form stars and clusters of low and intermediate mass, 
just as Ophiuchus and Perseus do. 

Still, many IRDCs will turn towards MSF. Interest- 
ingly, Table Q] suggests that most of the mass located 
in IRDCs is in clouds that will form massive stars. 
For example, th e 250 most compact clouds from the 
iPeretto fc Fullerl (|2009f ) sample (identified assuming a 
common distance) contain more than 50% of the area- 
integrated column density of all IRDCs. This suggests 
that they also contain a major fraction of the mass seen 
in IRDCs. If this reaso ning is correct, j ust few 10 2 
IRDCs (and not all - 10 4 : IRathborne et al.ll2006D might 
contain most of the Galaxy's star-forming gas. Given 
the uncertain na t ure of the properties derived from the 
IPeretto fc Fullerl (|2009( ) data (Section [373]), this conclu- 
sion is far from certain, though. 

5. CONCLUSIONS 

This letter studies whether Infrared Dark Clouds 
(IRDCs) are able to form massive stars. Our main con- 
clusions are as follows. 

• Observations of regions with and without massive 
star formation (MSF) suggest that the condition 
m(r) < 870 M© (r/pc) 133 (Eq. QJ approximates a 
threshold for MSF (Section El]). MSF clouds differ 
from those obeying Eq. (fTJ in mass-size slope or 
intercept (Fig. [TJ Section 14.11) . 

• Many IRDCs (Section I3.2[) . if not most (Section 
13. 3j) . do not exceed Eq. ([T]). Without significant 
further evolution, such clouds are unlikely candi- 
dates for MSF, but they might well form stars and 
clusters of up to intermediate mass (like Perseus 
and Ophiuchus). Very opaque IRDCs of large an- 
gular size constitute good MSF candidates. 

• Provided extinction-based masses can be trusted, 
just few 10 2 IRDCs might contain a major frac- 
tion of the Galaxy's star- forming gas fSection l4.3l) . 
These IRDCs would be dense and massive enough 
to host MSF. 



We arc indebted to a careful referee, who significantly 
helped to improve the text. This research was sup- 
ported by an appointment of Jens Kauffmann to the 
NASA Postdoctoral Program at the Jet Propulsion Lab- 
oratory, administered by Oak Ridge Associated Univer- 
sities through a contract with NASA. 
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